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A modified form of non-locally corrected theory of gravity is investigated in the context of 
cosmology and the Newtonian limit. This form of non-local correction to classic Einstein- 
Hilbert action can be locally represented by a triple-scalar-tensor theory in which one 

of the scalar degree of freedom is phantom-like and the other two are quintessence-like. 
We show that there exists a stable de Sitter solution for the cosmological dynamics if a 
suitable form of potential function V{4>) (or equivalently, f{R)) is selected. However, no 
matter what a potential function is selected, there is always an early time rcpeller solution 
corresponding to a radiation dominated universe. Besides, the equations for linear scalar 
perturbations are presented and it is shown that the form of potential function V{<p) is 
stringently constrained by the Solar system test, although the post-Newtonian parameter 
7 is not directly ailected by this function. 

PACS Nos.:98.80.-k, 95.36.-fx, 04.50.-h 
1. Introduction 

Strong evidences from the current cosmological observations converge upon the fact 
that the universe is spatially flat and there exists exotic component, dubbed dark 
energy, which drives the speed-up expansion of the universe. Many scenarios have 
been proposed to explain the acceleration in the framework of general relativity 
(GR). The preferred and simplest candidate for dark energy is the Einstein's cos- 
mological constant which can fit the observations well. However, it suffers from the 
so-called fine-timing problem and coincidcince problem. Meanwhile, the observations 
are not yet able to confirm that dark energy is indeed a constant. Actually, many 
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dynamical models of dark energy have been studied extensively during the past over 
ten years, for a comprehensive review, see for example HH 

As is pointed out, the dark energy problem may be essentially an issue of quan- 
tum gravity Although a complete theory of quantum gravity has not been estab- 
lished, some valuable ideas are thought to be the features of the theory of quantum 
gravity, such as non-locality. Non-locality can naturally arise in string theory and 
the non-local corrections to the effective action of gravity can be generated from 
quantum loop effects. It is also believed that the black hole information paradox 
may be resolved by the effect of non- locality l^l. Recently, a new non-local corrected 
theory of gravity was suggested in the form 

S^j^Jd^xR{l + f{a-^R)), (1) 

where D~^R denotes the inverse d'Alembertian acting on the scalar curvatur^. In 
the limit of / = 0, it goes back to general relativity. It was found that this kind of 
non-local gravity may lead to a unified description of the early-time universe and the 
present accelerating universd333 and there is no apparent conflicts with observations 
in the solar system test^. Non-local gravity has been extensively investigated in the 
literature I«|iu|ii|i^|i^|i4|i5|ie|i7|i8| 

In this work, we shall suggest a modified non-local term to the effective action 
of gravity and investigate its implications in cosmology and the solar system tests. 
Although, just as that in action ([l]), this form of non-local corrections to gravity 
has also no direct support from fundamental theories, we wish it appears to pro- 
vide a small window through which one can look into phcnomenological aspects of 
fundamental theories. 

This paper is organized as follows: we first present the model in general form 
and a triple-scalar tensor representation is given in Sec. [2j In Sec. [3j we investigate 
the cosmological dynamics of the theory. In Sec. |4] the scalar mode of cosmological 
perturbation to the first order and its Newtionian limit are considered. Finally, in 
section [5] we summarize our results and give some discussions. 



2. Non-local f{R) gravity and triple-scalar-tensor theory 

We here consider the following simple modified action of the non-locally corrected 
gravity 

5 = J d^xV^!^^[R{l + n-'f{R))]+C,nY (2) 

where g is the determinant of metric tensor g^^, k = SnG and / is a generic 
function which has the same dimension with Ricci scalar R. Here, \3~^ is the inverse 
d'Almbertian operator, which is in fact an integral operator, therefore, the above 
action is non-local. The equations of motion derived directly from this action will 
be rather difficult to handle because of the non-locality. However, by introducing 
two Lagrange multipliers 77 and ^ and replacing the inverse d'Alembertian and Ricci 
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scalar as two auxiliary scalar fields ijj and if respectively, the above action can be 
rewritten in the following form, 

-C(nV'-/M)} + ^™}- (3) 

It can be checked, by varying action ([s]) with respect to 77, ^ and respectively, that 
R = (fi, Dip — f{(p) as we expect and rj = 1 + ip + £,f'{(p). Integrating partially and 
neglecting the boundary terms, then the double derivative is eliminated as usual, 
we obtain a kind of triple-scalar-tensor theory 

5 = 1 d'xV^ |^[i?(l + ^ + e/'M)-^a'M 

+d^^dxiP + ^f{iP)] (4) 

It is noted that this kinetic coupling can not eliminated by a field redefinition. To 
be convenient, defining (j> = f'{(p) and taking ip as a. function of (j), then introducing 
a new function 1^(0) — 'p{4>)<p ~ fi'^{4'))i we can rewritten action Q as 

I d^x^l^^[R{l + iP + ^cP) 

-^V{^)+d^^dxi^]+C^}. (5) 

Contrary to the case for the action ([T]) , in which the local equations of motion can 
be derived by varying with respect to the metric and two auxiliary fields, there are 
three auxiliary local fields in our model ([s]). 

By varying action([5]) with respect to the metric g^^, we obtain the following 
modified Einstein equations 

+ (5^,n - V^V,)F - AtT^„ (6) 

where we have defined F — 1 + ip + S^cf) and G^^ is the Einstein tensor defined as 
usual and T^i, denotes the stress tensor of matter. Furthermore, the variation with 
respect to ip, and (p gives the following equations of motion 

□^ = y'(0)0-F(0), (7) 

□C = V'icb) (8) 

and S,{R — V'{(p)) — 0. When ^ is not zero, we have 

R = (9) 

where the prime denotes the derivative with respect to (p. Because matter is mini- 
mally coupled with other fields, the matter energy-momentum tensor is conserved, 
'^^TJJl = 0. Note that this conservation law is not an independent condition. 
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3. Cosmological dynamics 

In this section, we shall investigate the cosmological dynamics of the non-locally 
corrected theory presented in the above section. In a flat 4-dimensional Fricdmann- 
Robertson-Walker(FRW) space-time, 

ds^ = -dt^ + a{tf [dr"^ + r^de"^ + sin^ ed(j)'^] , (10) 

we obtain the Friedmann equations 

ZFH-'-^{iV{(l>)-ii>) + iFH^np^, (11) 

F + 2HF + {2H + 3H^)F - ^(^1^(0) + C^) = -upm (12) 

where H — d/a is the Hubble parameter, the overdots denote the derivatives with 
respect to cosmic time t, p,„ and P„j are energy density and pressure of matter, 
respectively. 

The two scalar field equations become 

ij + 3Hi^ + V' {((>)((> -V{(f>)^0, (13) 

i+3Hi + V'{(P) =0. (14) 

Besides, the usual matter conservation law reads 

p;„ + 3i/(l + w)p^ = 0, (15) 

where we have assumed an equation of state for the matter p„j — wpm- It is worth 
noting that, due to the Bianchi identity, only four equations among Eqs.( [Tlp5l ) are 
independent. 

To rewrite the equation of motions as a dynamical system, we may define new 
dimensionless variables 

y(</)) l/'(0) 

as the functions of — In a. Then, the cosmological dynamics can be represented 
by the following equations 

x\=X2, (16) 



2x3x7 + X4,xe + xyxs(1 + xi + X5 — 3 — ) 

2^3 



+ (^2^72^8 - ^)X2 (3W - l)fim, 

6 



(17) 
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X4, 



{-X7Xs - 1)X4 + XjXs, 
D 



X6, 



^6 = il^rxs - l)a;6 + (i^^ + l)x7, 

6 X3 



ilx,xs + A)x,-'-^ix,x,-x,X4), 



x'g ^ -(r - 1)-^{X2X3 - X1X4), 



(18) 
(19) 
(20) 
(21) 

(22) 
(23) 



-x-jx^ + 4 — 3(1 + w) 



(24) 



where the primes denote the derivative with respect to N and the parameter T is 
defined to be F = VV" /V'"^ . The Friedman equation is now a constraint equation 



1 1 

i^m = 1 + 2^1 + 2^2 + + + -p.XiX(i - -X3X7. 

6 6 



(25) 



For an exponential potential 1^(0) = Vbe"*", F = 1. Then, from Eq.(23), wc know 
xs = —a is a constant. One the other hand, it is clear that 2:7 = is a solution of 



Eq.(22|. Therefore, we can easily obtain the following solution 



X2 — C2e 



-N 



X6 = Cge 
X5 = C5- Cee"^, 

-N 

X4 = 046 , 
X3 = C3-C4e~^, 
(3w - l)^m ^(l-3m)jV _ 1 

3w-2 3 



-2N 



Xi — Ci — 026 — -e^ + -04066 , 

6W — 2 6 

where C.j(i = 1..6) and r2^„ are constants. The Friedman equation requires that 
1 + Ci + C5 = 0. From xj — 0, we know in this solution the scalar curvature R — 0, 
hence, a = aoVl — 2iJoio + '^H^t, where we have set H{ti^) = Hq and a(io) = Uq. 
Clearly, when the equation of state w > 1/3, the solution is stable and the fixed 
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point {xi — Ci,X2 — 0,X3 — C3,X4 — 0,X5 — C'c, — ~Ci — 1,xg = 0,X7 = 0) is 
an atrractor. In figure^ we show the evolution of effect equation of state w^ff = 
— — 1 in the exponential model for different initial conditions. It is found that 
there exist a stable attractor where We// = 1/3. We note that in this model, the 
acceleration of the universe depend heavily on the ininital conditions. This solution 
may be a good representation of the early universe when a radiation dominated 
phase follows naturally from an inflationary expansion. Also, it is interesting to find 
that there is a case with realizing a crossing of the phantom divide from the non- 
phantom (quintessence) phase to the phantom one for some time and then going 
back to the non-phantom phase at late time. 




Let us now consider whether there is a de Sitter solution of the theory in which 
the Hubble parameter is a constant, H{t) — Hq and the scale factor a ~ aoe^"*. 
We assume that ^ is not zero, so V'{(l)) in the equations is just the scalar curvature 
R, which is also a constant in this case, R — ISiJg. Therefore, from Eq.(14|, the 



4, 2011 1:29 WSPC/INSTRUCTION FILE mainV2 



A modified model of non-locally corrected gravity: cosmology and Newtonian limit 7 

solution for ^ can be straightforward obtained, 

^ = C2^AHot-^e-^"°K (26) 

For simplicity, we let the intergal constants Ci and C'2 vanish. This means 2:3 = 
—AN. We also assume the potential function V{(l)) is given by a linear function 
V{(j)) = Vq4), where the constant Vq should be equal to 12Hq due to the relation 
([9]). Therefore, from the definition of x-j and xs, we have xtXs = —12. Applying 
this relation and assuming X2 — 0, we obtain that a;i = A is a constant, X4 — —4, 
X5 = A-l-^Be-^^, xe = Be-^^ , xr = -3A/7V, xs = 4iV/A, = 17°„e-3(i+"')^, 
where A, B and 51^^ are three constants. Thus, using cosmic time t, we can represent 
the solution as = -41^, F = 2A - fe-^^"* and V = ^ - 1 - f e-^^"*. Note 
that if the potential y is a linear function of 0, from its definition, f{R) should 
be also a linear function, that is, the non-locally corrected term in the action reads 
R\Zr^[R — Vq). In Fig. [2[ we plot the trajectories of the evolution of We// under 
different initial conditions. It is shown that Weff trends to the limit —1 in the 
future regardless of what the initial value takes, which may means that the de Sitter 
solution above is a stable attractor solution. Meanwhile, in the past, there may exists 
an rcpellcr solution representing an effective radiation dominated universe in which 

Weff = 1/3. 




Fig. 2. The evolution of the effect equation of state w^jf with respect to Ina in the linear model 
y(0) = Vo<l> for different initial conditions. 
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4. Scalar perturbations and Newtonian limit 

In this section, we show cosmological perturbation equations in our theory. Clearly, 
the evolution of scalar perturbations are much more complicated than those of the 
vector and tensor perturbations because they are heavily coupled to scalar fields 
and matter. Here we only consider scalar perturbations. As usual, once the scalar 
modes of perturbations are introduced, the perturbed metric can be written in the 
Newtonian gauge to the first order as 



(1 + 2^{x, t))dt^ + a(tf5^j[\ + 2$(f, t))dx'dx\ 



(27) 



where 5* and $ are two scalar potentials. When both of the two potentials vanish, the 
FRW metric of the background universe is recovered. Meanwhile, the three scalar 
fields are perturbed to be ip{x,t) — '4){t) + S'ilj{x,t), S,{x,t) = ^(i) + S£,{x,t) and 
^(x, t) = (pit) + 54>{x, t) respectively. Similarly, as for matter, the energy density 
and pressure are perturbed to be pm{x,t) — pm(t) + 5pm{x,t) and Pm(x,t) = 
Pm {t) + 5Pm (x, t) . Inserting these expressions into Einstein equation ([6| , up to the 
first order, we obtain the following equations for the perturbations 



2F 
1 

+ 2 



^3H(2*F - SF) - 3<i>F + a-^V'SF = -nSp, 



(28) 



and 



+{SF)^o,-F^^,-H{6F)_,^0 



(29) 



S% { {2H + 3ij2)(2F* -SF) + - [V{<j))S^ + ^RSc/)] - 2F$ 

2 

+ (* - 3$)(2iJF + F) + 2*(F + 2HF) - {SF + 2H6F) 
+ a-^ [F(V2$ + V^*) + \/^5F] } 



F{^ 



a-^SF, 



^S'jSPra, 



(30) 



where SF — 6'ip + (j)d^ + ^S(j3. Similarly, from equations of motion for the scalar fields, 
we have 



{6^ + 3HS^) + ^(* - 3$) - 6 $ - H{if - 4$) 



2* 



a-^ [\7^5£. + 2(V2* + 2\/^^)] = 0, 



(31) 
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and 



- {S'i^ + SHS't/j) + - 3$) - 60 $ - - 4$) 
+ 2* + SHiP) + 6(j){H + 2H^) 

+ a-^ [V^S^j + 2(/)(V2* + 2V2$)] = 0. (32) 

In the Newtonian Hmit of the cosmological perturbations, we may omit the 
cosmological expansion and set scale factor a{t) = 1. Moreover, the derivatives with 
respect to time, relative to those with respect to space, are also negUgibly small. 
Therefore, in the vacuum outside a mass where the source term can be neglected, 
the shear constraint (|30|) gives 



o 



which shows that the gravitational potentials are not equal, 

SF 



* + $ = - 



F 



Meanwhile, from the perturbation equations (31) and (32), we have 

V^S^ + 2{VH + 2V2$) = 0, 
V^J?/' + 2(t){V'^^ + 2V2$) = 0. 



(33) 

(34) 
(35) 



2 



(36) 



Solving the above two equations and using the condition that the usual 
Schwarzschild solution is recovered in the limit of general relativity, we obtain that 

The perturbation 5(j) may be neglected because field <}> is not dynamical. Therefore, 
applying the relation 5F = Sij: + (f)S^ to eliminating 6F, 5^ and Scf), we get the 
post-Newtonian parameter 

$ _ - 40 
^ ~ F-Sch' 



7 



(37) 



which is stringently constrained to its value in general relativity j — 1 from the 
solar system tests such as the Cassini spacecraft It should be noted that 7 is 
not directly affected by the potential function 1^(0) (or equivalently, f{Rj), the 
constraints is still imposed because the value of F and are heavily dependent on 
y(0). 



5. Conclusions 

In this work, we investigated a new form of non-locally modified theory of gravity 
with action (|2|. Just as action ([T]) can be recast as a local biscalar-tensor theory, 
this non-local theory can be locally represented by a triple-scalar-tensor theory. In 
the context of cosmology, we show that there exists a stable de Sitter solution for 
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the cosmological dynamics if a suitable form of potential function V{4>) (or equiva- 
lently, f{R)) is selected. But, no matter what a potential function is selected, there 
is always an early time repeller solution corresponding to a radiation dominated 
universe. However, it should be stressed that a long enough period of matter dom- 
ination, in which structure is supposed to be formed, is crucial for a successful 
cosmology and a lack of a matter dominated era could make an accelerating model 
impracticable. Therefore, it may be an alternative to dark energy for the late-time 
acceleration, only if a matter dominated era can appear naturally between the ra- 
diation dominated era and the late-time accelerating phase. On the other hand, we 
presented a set of equations for scalar linear perturbations and considered the New- 
tonian limit for static spherically symmetric solutions. It is found that the form of 
potential function V{(t>) is stringently constrained by the solar system test, although 
the post-Newtonian parameter 7 is not directly affected by this function. 

Just as in the framework of most higher derivative gravity theories it has 
been shown that ghosts will appear in most of the non-local gravity models. In 
Ref.'^, it is also found that the introduction of a f{R) term into the non-local 
action ([l]) will not contribute to the solution of this problem. Unfortunately, in the 
non-local theory investigated here, this problem is still unsolved. One of the three 
introduced scalar degree of freedom is phantom-like (or a ghost) and the other 
two are quintessence-like. See Appendix A for the detail discussions. Introducing 
a new non-minimally coupled scalar field into the action of gravity may alleviate 
this problem ^221. Although the theory considered here may have some problems in 
some aspects and there is no direct support from fundamental theories, it is hard 
to deny that it appears to provide a small window through which one can look into 
phenomenological aspects of fundamental theories. 

As is pointed out in Ref. ^^^l^ replacing the operator in action ^ will intro- 
duce spurious degrees of freedom, because even when f{D^^R) is constant, there 
exist nontrivial homogeneous and isotropic vacuum solutions of the Einstein solu- 
tion in the scalar-tensor representation of the theory. However, we find that the two 
cosmological solutions considered in the work are not artifacts of the scalar-tensor 
theory ([5|, but should be two real solutions of the original theory ([2|, although we 
are now not sure that whether the field content of the original theory and the triple- 
scalar-tensor representation are the same. This may deserve thorough investigation 
in the future. 

Appendix A. Einstein frame and ghost problem 

Let us perform a conformal transformation to the Einstein frame 



(A.l) 



hence 



R = n-^ [R^ - 6(n \nn + g^'^V^ In 17V In n)] , 



(A.2) 



4, 2011 1:29 WSPC/INSTRUCTION FILE mainV2 



A modified model of non-locally corrected gravity: cosmology and Newtonian limit 11 

with = (1 + ■;/; + ^0)^^ = e*. Therefore, for the vacuum case, action ^ can be 
rewritten as 

S = [ d^x^ [r - ^V^$Va$ - e*V^eVA$ 
J 2n I 2 

- e^^cf^^^Vxi - e2*ev^CVA0 - e^^^V{4>) ] (A.3) 

Defining M as the coefficient matrix of the kinetic terms 

/ I ie* \ 
M = Ue* e2*0 i^e"* (A.4) 

and A as its eigenvalue. Since M is real 3x3 symmetric matrix, there must exist 
three real eigenvalues for the kinetic terms, which are determined by the condition 
det(M — IX) — 0. In order to avoid the ghost degree of freedom, all the three values 
of A must be positive. However, the product of the three values is determined by 
det M — — fe"'*^^, which is always negative and the sum of them is | + e^*(/) which 
is positive. This means that one of the three eigenvalues of the coefficient matrix 
M is negative and the other two are positive. 
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